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ABSTRACT 



The Fermi Large Area Telescope (LAT) data have confirmed the pulsed emis- 
sion from all six high-confidence gamma-ray pulsars previously known from the 
EGRET observations. We report results obtained from the analysis of 13 months 
of LAT data for three of these pulsars (PSR J1057-5226, PSR J1709-4429, and 
PSR J1952+3252) each of which had some unique feature among the EGRET 
pulsars. The excellent sensitivity of LAT allows more detailed analysis of the 
evolution of the pulse profile with energy and also of the variation of the spectral 
shape with phase. We measure the cutoff energy of the pulsed emission from 
these pulsars for the first time and provide a more complete picture of the emis- 
sion mechanism. The results confirm some, but not all, of the features seen in 
the EGRET data. 



Subject headings: gamma rays: stars — pulsars: individual (PSR J1057— 5226, 
PSR J1709-4429, and PSR J1952+3252) 
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Introduction 



Gamma-ray pulsars offer valuable probes of rotating neutron stars. The gamma-ray 
emission is directly related to the primary particle acceleration processes in the pulsar mag- 
netosphere. This high-energy (E^ >100 MeV) radiation represents a significant fraction of 
the spin-down luminosity of p ulsars in s ome c ases. For a recent summary of implications of 
gamma-ray pulsar studies, see lHardingi ( 120081 ). 



Until the launch of the Fermi Gamma-ray Space Telescope (Fermi) on 2008 June 11, 
the most complete results on the highest-energy emission from pulsars came from the En- 
ergetic Gamma Ray Experiment Telescope (EGRET) on the Co mpton Gamma R ay Ob- 
servatory, which found high- confidence detections of six pulsars ( [Thompson! 120081 ) . The 
Larg e Area Telescope (L AT) on Fermi has detected a large population of gamma-ray pul- 
sars (lAbdo et al.l l2010el ). including all six of the EGRET pulsars. Detailed st udies of the 



three brighte st of th ese have been repo rted by the LAT Collaboration: Ve la (lAbdo et al. 
2009al . boioj l. Crab (Ubdo et al.l boiOch . and Geminga (Ubdo et al.l hoiOdh . The present 



paper is a study of the remaining three EGRET high- confidence pulsars, PSR J1057— 5226 
(B1055-52), PSR J1709-4429 (B1706-44), and PSR J1952+3252 (B1951+32). The key 
measured and derived parameters of these pulsars are listed in Table [lj 



Pulsed emission from PSR J1057— 5226 was first detected in radio by lVaughan fc Large 

(119721 ) with a period of 197.11 ms. It belongs to th e small class of r adio pulsars with 
a strong interpulse midway between the main pulses (IKeith et al.ll2010l ). No evidence of 
pulsa tion has been observed in optical, but it has been detected by the Ein stein Observa- 



1993) and it has been intensively studied with Chandra ( iTeter et al.l 120011 ) . iBecker et al. 



tory ( jCheng fc Helfandlll983l ) and by EXOSA T (IBrinkmann fc Ogelmanlll9871) as a soft X-ray 
sourc e. In 1993, pulsed emission in the X-ray band was detec ted by ■RQ&4r(l0gelman fc Finlev 



19991 ) described the unsuccessful search for a pulsar wind nebula around this source. 



The geometry of this pulsar has been studied in detail by IWeltevrede &: Wrightl ( 120091 ). 
based on Parkes radio observations. Fitti ng the phase sweep of the pol arization position 
angle to the rotating vector model (RVM) (IRadhakrishnan fc Cookelll969l ). they found that 
the pulsar's magnetic axis is inclined at an angle of ~ 75 to its rotation axis and that the 
radio interpulse arises from emission formed on open field lines close to the magnetic axis 
that do not pass through the magnetosphere's null (zero-charge) surface. The radio main 
pulse emission must originate from field lines lying well outside the polar cap (PC) boundary, 
beyond the null surface and farther away from the magnetic axis than those of the outer gap 
(OG) region where the gamma-ray peak is generated. 



PSR J1709— 4429 was first detected as an unidentified gamma-ray source by COS- 



-6 - 



B (jSwanenburg et al.lll98ll ). Its identification as a pulsar was made abou t 10 years later 



with the discovery of pulsed emission from this source in radio wavelengths (IJohnston et al. 



19921 ) and wi th the detection of the gamma-ray pulsations by EGRET using the radio timing 
informatio n (IThompson et al. I ll996h . The pulsed emission in X-rays was found in 2002 with 



Chandra (Gotthelf et al. 



in radio (IFrail et al. 



1994 



2002h. A diffuse pu lsar wind ne bula was seen around the pulsar 



Giacani et al.ll200ll ) and X-rays ( iFinley et al.lll998r iRomani et al. 



20051 ) with an extension of 3' and 110", respectively. A possible association with a faint su- 



pernova remnant (SN R) G342. 1—2.3, imaged as an arc-like radio structure (IMcAdam et al. 



19931 ; IFrail et al.lll994j ). was also suggested, although the implication of a high proper motion 
velocity due to the off-center location of t he pulsar in this stru cture conflicts with the mea- 
sured scintillation velocity of the pulsar ( INicastro et al.lll996l ). No point-like emission was 
detected in the very-high-energy (VHE) band from the HE SS observations at the location of 
the pulsar (lAharonian et al.l 120051 ; lYoshikoshi et al.l 120091 ); however, extended emission was 
detected from a location offset by ~0.35° from the pulsar, consistent with the SNR radio 



location with an extension of 0.3° (IHoppe et al.ll2009l ). 



The 39.5 ms ra dio pulsar PSR J1952+3252 in the CTB80 SNR was disc overed by 
Kulkarni et al.l ( 119881 ). Proper motion measurements of the pulsar using the VLA0 Jzdg eretaJj 
20081 ) leads to a Local Standard Corrected transverse velocity of 274 km s _1 and a kinetic age 



of 51 kyr assuming the pulsar was born in the geometric center of CTB 80 and a distance of 2 
kpc. Pulsations in th e X-ray domain have been reported from EXOSAT (jOgelman fc Buccheri 
19871 ) and ROSAT (IBecker fc Triimperl Il996l ) with low signifi cance. In the GeV range, a 



pulsed emission showing two peaks was observed by EGRET (IRamanamurthv et a 



No emission has been detected at higher energy ( Zweerink et al.ll2009l ; lAlbert et al. 



1995). 



20071 ). 



Each of these pulsars had some unique features among the EGRET pulsars that merit 
investigation with the greater sensitivity, better resolution, and broader energy range of 
the LAT. PSR J1057— 5226 in the EGRET data appeared to ha ve the highest efficienc y 
for conversion of spin-down luminosity into gamma-ray energy (IThompson et al.l 119991 ). 
PSR J1709— 4429 had a measured EGRET energy sp ectrum that was well de scribed by a bro- 
ken power law rather than a sharp high-energy cutoff (IThompson et al. Il996h . PSRJ1952+3252 
showed an energy spectrum in the EGRET data with no evidence of a cutoff or spectral 
break out to the 30 GeV limit of the EGRE T energy range, although EGRET found only 
two p hotons above 10 GeV from this pulsar (IRamanamurthv et al.lll995l ; IThompson et al. 



2005) 



x The VLA is operated by the National Radio Astronomy Observatory, which is a facility of the National 
Science Foundation, operated under cooperative agreement by Associated Universities, Inc. 
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In addition to updating conclusions from past observations, the detailed LAT observa- 
tions of these pulsars have intrinsic interest. PSR J1057— 5226 and PSR J1709— 4429 are 
among the minority of gamma-ray pulsars having double peaks in their light curves with a 
separation less than 0.3 of the phase. PSR J1952+3252 is a short-period pulsar and lies in 
an SNR/pulsar wind nebula complex. 



All three of these pulsars are included in the first Fermi LAT pulsar catalog ( lAbdo et al. 



2010eJ). The present work, in addition to including more data than the catalog paper, also 



extends the analysis in several ways: (1) a more detailed analysis of the pulse profiles as a 
function of energy; (2) phase-resolved spectra; (3) a search for off-pulse emission (such as a 
pulsar wind nebula); and (4) a more detailed discussion of models for beaming and efficiency 
calculations. 



Observations 



2.1. Gamma-ray Observations 



The LAT aboard Fermi is an electron-positron pair conversion telescope sensitive to 
gamma rays with energies in the range from 0.02 to greater than 300 GeV. The LAT is made 
of a high-resolution silicon micro-strip tracker, a Csl hodoscopic electromagnetic calorimeter, 
and a segmented pla stic scintillator detector to identify the background of charged particles 
( lAtwood et al.ll2009[ ). Compared with its predecessor EGRET, the LAT has a larger effective 
area (~8000 cm 2 on-axis, >lGeV) and improved angular resolution (# 68 ~ 0.6° at 1 GeV 
for events in the front section of the tracker). The large field of view (~2.4 sr) allows the 
LAT to observe the full sky in survey mode every 3 hr. The LAT timing is derived from a 
GPS clock on the spacecraft, and gamma rays a re hardware time-stamped to an accuracy 
significantly better than 1 [is ( lAbdo et al.ll2009bl). The LA T software tools for pulsars have 
been shown to be accurate to a few /xs (ISmith et al.l 120081 ) . 



2.2. Radio Observations 



These three pulsars have been monitored by several observatories as part of the pulsar 
timing c ampaign for Fermi , monitoring the E > 10 34 erg s _1 pulsars at radio and X-ray wave- 
lengths jSmith et al.l l2Q08h . PSR J1057-5226 and PSR J1709-4429 are observed monthly 
with the 64 m Parkes radio telescope in Australia. Typical observations last for 2 minutes at a 
frequency of 1.4 GHz (and occasional obs ervations at 0.7 and 3.1 G Hz). Full details of the ob- 
serving and data analysis can be found in lWeltevrede et al.l (120101 ) . Both pulsars have a high 
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->35 



degree of linear polarization in the rad io, as is often the case for pulsars with E > 10 — 10 
ergs -1 (jWeltevrede fc Johnston! 120081 ). The timing solution for P SR J1952+3252 uses radio 
observations made at the Nangay telescope (jCognard et al.ll2009l ) at 1.4 GHz. 



3. Timing Analysis 



The radio timing solutions for these three pulsars have been derived from multiple 
times of arrivals (TOAs) using TEMP02 (IHobbs et al.ll2006l ). The TOAs were fit to give 
each pulsar's spin frequency and frequency derivative. The data are whitened using the 
fitwaves algorithm within TEMP02 in order to take into account timing noise. The number 
of TOAs used in building the timing solutions, the validity period of the solution, and the 
resulting timing rms for each pulsar are listed in Table [2j The complete timing solutions 
will be made available at the Fermi Science Support Center (FSSC) Web site§. When more 
than one timing solution is available, we have used the one with the lower rms value in order 
to construct a pulse profile as fine as possible and not to be limited by the accuracy of the 
timing solution. 

In the 13 month interval defined in Section HI a glitch, a sudden change in the rotation 
speed of the pulsar, was detected in the analysis of the LA T data for PSR J170 9— 4429. 
While a glitch has been detected from this pulsar previously ( IJohnston et al.lll995f ). this is 
one o f a few detected and characterized by a gamma-ray telescope using gamma-ray data 
only (j Jackson et al.l 120021 ; ISaz Parkinson! 120091 ) . The time of the glitch has been narrowed 
down to a ~22 hr window (MET@ 240395907-240472803) around 2008 August 14 with the 
LAT data. The parameters of the glitch were measure d from the radio data as Au/u = 
2.7497 ± 0.0001 x 10~ 6 and Az>/z> = 4.95 ± 0.01 x 10 -3 JWeltevrede et al-lboid ). 



PSR J1709— 4429 was also timed with the LAT data alone obtained over 13 months of 
LAT observations. The timing solution derived from the LAT observations fits the frequency, 
frequency derivative, and glitch parameters for this pulsar as well as six WAVE terms to 
whiten th e stron g timing noise of this pulsar, using the analysis techniques described in 
Ray et al.l (120101 ). We used the radio timing solution with the LAT data to define the radio- 
gamma offset in the pulse profile constructed with the LAT timing solution. 



2 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/ephems/ 

3 MET is the Mission Elapsed Time in seconds from 00:00 UTC on 2001 January 1. 
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4. Fermi LAT Analysis and Results 

Fermi LAT data were analyzed using the standard Science Tools (ST) in its v9rl5p2 
version available at the FSScfl Only "Diffuse" class events, which have the tightest back- 
ground rejection, were selected for the analysis. In addition, we excluded those photons 
coming from zenith angles > 105°, where the gamma-rays from Earth's limb produce exces- 
sive background contamination. We also have excluded time intervals when the region of 
interest (ROI) intersects Earth's limb. 

The data set spans 13 months from the start of the sky survey observations, 2008 August 
4 (MJD 54682.66) to 2009 August 28 (MJD 55071.94). We discarded the observations 
performed during the Launch and Early Operations (L&EO) period in order to avoid the 
effects of the instrument r esponse function (1 RF) changes due to various configuration tests 



during the L&EO period (lAbdo et al.ll2009bl ). The last date of all data sets was determined 
by the latest available ephemerides provided by the radio telescopes for each of the pulsars 
reported in this paper. Thus, the data set for each pulsar covers a slightly different time 
range. 



4.1. Pulse Profiles 

Like all pair-production detectors, the LAT's angular resolution is dominated by multiple 
scattering at low photon energies. It is therefore important to use an energy-dependent radius 
similar to the LAT point spread function (PSF) in order to include only the events that can be 
strongly associated with that source. The accurate parameterization of the LAT PSF to be 
used for science analysis is described by the IRFs. A simplified, acceptance-averaged approxi- 
mation of the PSF with a 68% containment angle is given by < 9 68 (E) >= (0.8°) x (E Ge y)~ S . 
Accordingly, for the pulse profile analysis, we selected only those events within an energy- 
dependent radius of 9 < Max(Min(i? max , # 68 ), 0.35°) around each pulsar: th e minimum 



value of 0.35° was set in order to be sure to keep all high-energy photons ( lAbdo et al. 



2010ef ). while a maximum radius, -R max , was introduced to reduce the background contam- 



ination at low energies. Th e value of i?, max was determined such that it maximizes the 



if -test (Ide Jager et al.lll989f ) performed on the light curve built with only low-energy pho- 
tons (0.1 GeV < E < 0.3 GeV): this analysis found i? max of 2.0° for PSR J1057-5226 and 
PSR J1952+3252 and 2.4° for PSR J1709-4429. Note that this ROI cut was only used to 
select data for obtaining the pulse profiles, but was not applied to the data sample used in 



4 http://fermi. gsfc.nasa.gov/ssc/data/analysis/software/ 
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the spectral fits. 

Finally, photons with energies above 100 MeV were used for this analysis. The energy 
cut is necessary because all three pulsars are close to the Galactic plane, and due to the poor 
angular resolution below 100 MeV the photons from the strong Galactic diffuse emission hide 
any periodicity. 

In the following subsections, we report the pulse profile characteristics for each pulsar. 
We give the total number of pulsed and background photons obtained with the selection 
above. For each pulsar, we show a full-band pulse profile and its evolution with energy, 
constructing pulse profiles for four different energy bands of 0.1— 0.3GeV, 0.3— l.OGeV, 
1.0-3.0 GeV, and >3.0GeV for PSR J1057-5226 and PSR J1952+3252; and for five en- 
ergy bands of 0.1-0.3 GeV, 0.3-1.0 GeV, 1.0-3.0 GeV, 3.0-10.0 GeV, and >10.0GeV for 
PSR J1709— 4429. These pulse profiles have been built with a fixed bin width of 0.01 in 
phase except for the highest energy band for which we have used a bin width of 0.02 be- 
cause of limited statistics. In these gamma-ray profiles, ip = is defined as the maximum 
of the main radio peak for PSR J1709-4429 and PSR J1952+3252. The radio light curve 
for PSR J1057— 5226 exhibits multiple peaks and we have chosen to put the radio interpulse 
peak at phase zero. We hence shifted the gamma-ray profile accordingly. The pulse shapes 
in each band are parameterized by fitting with simple functions, and the evolution of these 
parameters with energy is discussed for each pulsar. 



Using the data selection and energy- dependent ROI described in Sections @] and I4.1[ 
we detected 6106 events with energies above 100 MeV in the region around the pulsar. 
Since there is no evidence of significant off-pulse emission (see Section 14 . 3 j) . we evaluated 
the number of background photons from the off-pulse phase region, and found 2834 ± 80 
background photons and 3272 ± 112 pulsed photons. 

The pulse profiles in different energy bands are shown in Figure [1] for two complete 
pulsar rotation cycles. The complex pulsed emission profile extends from 0.25 to 0.65 and 
changes significantly with energy. Although they do not represent a major portion of the 
emission, small leading (PI) and trailing (P2) peaks are always visible. We fitted only these 
two peaks with Gaussians over a constant background. The position of PI is 0.31 ± 0.01 
and its half-width is 0.04 ± 0.01 while the position of P2 is 0.59 ± 0.01 with a half width of 
0.03 ± 0.01. The location and half-widths of PI and P2 obtained from the fits in all energy 
bands are given in Table El These two peaks are separated by 0.28 ± 0.03 in phase. The 



4.1.1. PSR J 1057- 5226 
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region between these two peaks is complex. More photons will be needed to resolve the 
shape. The positions of the two peaks are constant with energy within the fit uncertainties. 
Figure [2] shows the P1/P2 ratio, seen to decrease with increasing energy. From Figured! 
P2 is 0.03 ± 0.01 after the peak of the radio main pulse and PI is 0.75 ± 0.03 after it. We 
find 333 photons with energies greater than 3 GeV, and the highest energy photon consistent 
with the pulsar position comes from the trailing peak, with a phase of 0.61 and an energy 
of 8.7 GeV. 



4.1.2. PSR J1709-U29 

Due to timing uncertainty around the glitch, we have chosen to exclude all the photons 
detected before the end of the glitch relaxation time on 2008 September 16 (MET 243216004). 
Using the data selection and energy-dependent ROI described previously, we detected 30448 
events with energies above 100 MeV in the region around PSR J1709— 4429. Failing to detect 
any significant off-pulse source in the phase region <p = 0.7 — 1.0 (Section l4.3p . we estimated 
the number of photons due to the background as 17004 ± 130 from this off-pulse phase range. 

Figure [3] shows the pulse profiles in separate energy bands for two complete pulsar 
rotation cycles. The full-energy-band pulse profile shows a dominant two-peak structure 
as opposed to the EGRET pulse profile which could not differentiate between two or three 
peaks. The first peak (PI) is located at 0.242 ± 0.002 in phase and the second peak (P2) 
is located at 0.492 ± 0.004 in phase. Both peaks are strongly asymmetric with sharp outer 
edges toward the off-pulse region with wider and structured inner edges toward the bridge 
region; hence, both peaks were characterized with asymmetric Lorentz ian function s , whic h 



have different widths for the leading and trailing edges, as described in lAbdo et al.l ( 12010fl ). 
The Lorentzian half-widths for the outer edges are 0.027 ± 0.002 and 0.032 ± 0.003, and the 
Lorentzian half-widths for the inner edges are 0.072 ±0.008 and 0.118 ±0.013 for PI and P2, 
respectively. The location and half- widths of PI and P2 obtained from the fits in all energy 
bands are given in Table |H 

The locations of the first and second peaks (PI and P2) stay constant over the energy 
bands, but PI seems to disappear in the highest energy band, > 10 GeV. The ratio of the 
first peak to the second peak slowly decreases with energy, as seen in Figure HI The pulse 
profile shows additional peak-like features between the peaks in the three highest energy 
bands. A third peak (P3) in the bridge region at the phase of 0.35 is seen above 3 GeV, 
similar to what is observed in the Vela pulsar, and P3 becomes as large as P2 above 10 
GeV. However, no phase shift of P3 with energy is detected for this pulsar with the current 
statistics. Above 10 GeV, 45 photons were detected from the direction of the pulsar in the 
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pulsed region, (p — 0.1 — 0.7. The highest energy photon has an energy of 47.7 GeV and 
is detected at phase 0.31. The approximate radio-Pi lag was found to be 0.24 using the 
method discussed in Section [3J 

4.1.3. PSR J 1952+ 3252 

Using the data selection and energy- dependent ROI described in Sections 0] and 14.11 we 
detected 10314 photons with energies above 100 MeV in the region around this pulsar. The 
number of background photons, estimated assuming no emission from the off-pulse region, 
ip = 0.8 — 1.0, as discussed in Section I4T31 is 8308 ± 159. 

Figure [5] shows the energy evolution of the pulse profile, which exhibits two clear peaks 
(further referred to as PI and P2). PI appears symmetric over all energy ranges whereas 
P2 is clearly asymmetric; accordingly PI is fitted with a symmetric Lorentzian and P2 with 
an asymmetric Lorentzian function over a constant background. The fitted functions are 
shown in Figure [5], and the location and half-widths of PI and P2 obtained from these fits 
in all energy bands are given in Table [51 The peak positions are independent of the photon 
energies and the FWHMs of the peaks decrease with increasing energy. Figure [6] plots the 
evolution of the P1/P2 ratio through the four energy bands defined previously, showing that 
PI fades with increasing energy whereas P2 persists. The highest energy photon in the whole 
sample is 25.7 GeV and belongs to P2 at tp = 0.58, and nine events with E > 10 GeV have 
been detected. 

PI lags the maximum of the 1.4 GHz radio peak shown in the bottom frame of Figure [5] 
by 0.154 ± 0.001. The separation of the gamma-ray peaks is 0.485 ± 0.017. 

4.2. Spectral Analysis 

For the spectral analysis of each pulsar, we used a data set that covers the same time 
period as the pulse profile studies and included all photons with energies in the range 
100 MeV — 100 GeV in the 10° region around each pulsar position. The "P6_V3" IRFs, 
which are a post-launch update to address gamma-ray detection inefficiencies correlated 
with the trigger rate, were used in the analysis. The systematic errors on the effective area 
are < 5% near 1 GeV, 10% below 0.1 GeV, and 20% over 10 GeV. In order to propagate the 
uncertainties on the effective area to the systematic errors on the three spectral parameters, 
the index, flux, and the energy cutoff, six "modified IRF" sets were used, each pair bracket- 
ing the nominal (P6_V3) effective area by a correction factor derived corresponding to each 
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spectral parameter. 

We measured a phase-averaged spectrum for each pulsar, and the results are reported in 
Section 14.2.11 All three pulsars were sufficiently bright that it was also possible to measure 
phase-resolved spectra for them. The phase-resolved spectra for these pulsars are reported 
in Section 14.2.21 Phase- averaged and phase-resolved spectra for the three pulsars are shown 
in Figures [7] and [8] for PSR J1057-5226, in Figures M and M for PSR J1709-4429, and in 
Figures HH and d for PSR J1952+3252. 



4-2.1. Phase-averaged Spectrum 



The emission spectrum for these pulsars in the energy range from 100 MeV to 100 
GeV was derived u sing the LAT tool qtlike, part of the ST package, based on a maximum- 
likelihood method ( iMattox et al.lll996l ). This tool fits a model representing the point sources 
and diffuse emission in the selected ROI to the data and finds the best-fit parameters to op- 
timize the likelihood function describing the data. The Galactic diffuse emission is modeled 
by the standard mapcube file gllJem_v2.fits. A tabulated model for the isotropic component, 
isotropic_iem_v02.txt, representing the extragalactic emission, as well as the residual instru- 
mental background, was used. Both diffuse models are publicly released by the FSSC. The 
absolute normalization was left as a free parameter for both mod els during the likel ihood 
fit. The nearby sources were taken from the First LAT Catalog (lAbdo et al.ll2010bl ). We 
included all sources within 17° of the pulsars. The spectra of all the nearby sources except 
pulsars were modeled with a power-law spectral shape. Pulsars were modeled as power laws 
with exponential cutoffs. The spectral parameters for sources within 10° of the pulsar were 
left free in the fit, and the spectral parameters for other sources were fixed to the spectral 
parameters obtained from the catalog run. 

In order to better evaluate the background in the region, we first fitted the diffuse 
background and nearby sources using the photons arriving in the off-pulse interval of the 
pulsars. The full-band phase-averaged spectrum for the pulsar of interest was then obtained 
using all the photons from the full phase range, and fixing the diffuse background and all 
nearby source spectral parameters to the best fit model obtained from the off-pulse analysis, 
properly rescaled to the full phase interval. At this step, the spectrum of the pulsar of 
interest was modeled with a power-law with exponential cutoff shape given by 
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where three parameters, the differential flux K (ph cm 2 s 1 MeV x ), the photon index T, 
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and the cutoff energy, Cutoff j were allowed to be free in the fit and 6 was fixed to be 1 (i.e., 
simple exponential). The fitted spectrum, with this assumed spectral model, for each pulsar 
is given in Table M and plotted as the solid curve in Figures EH and [IT] for these pulsars. 

The presence of cutoffs in the spectra of the pulsars under study was checked by fitting 
each spectrum with a power-law shape, and this model was rejected significantly for all three 
pulsars in favor of the cutoff model by a likelihood ratio test. A broken-power-law spectrum 
was the best-fit model for the PSR J1709— 4429 with EGRET data. This spectrum model 
was also tested for PSR J1709— 4429 with the LAT data, and it was rejected by 43.4a in 
favor of a simple exponential cutoff shape. 

Finally, the signal-to-noise ratios obtained from these pulsars were sufficiently high 
that we tested the validity of the simple-exponential-cutoff shape assumption by letting 
the 6 parameter free in the spectral fit model. We found that for PSR J1709— 4429 and 
PSR J1952+3252 the best-fit spectrum has b < 1, rejecting the 6=1 shape with 7.2cr for 
PSR J1709-4429 and with 2.9a for PSR J1952+3252. The fit for PSR J1057-5226 did not 
improve with a free b parameter and resulted in a value for b consistent with 1. The spectral 
parameters obtained from the fits with a free b parameter are also given in Table E] and the 
b < 1 fits are shown with the dashed curves in Figures [9] and [TTJ 

The flux points in Figures Ui and [11] were obtained with an analysis method that 
is independent of the model used to fit the sources in the region. We divided the data 
set into energy bins logarithmically spaced and in each of them we applied the gtlike tool, 
but this time, for all the point-like sources, we assumed a power law spectrum with a fixed 
photon index of 2, and a free flux parameter. We also modeled the Galactic diffuse emission 
and the isotropic component as before. In this way, we obtained the flux values for all the 
sources in each energy bin assuming that in the energy bin the spectral shape can be safely 
approximated with a power law. This method is more accurate as the energy bin size is 
reduced, but statistical considerations limit the minimum bin width. 

From the comparison of the spectral flux points and the full-band spectrum curve, it is 
clearly seen that while a simple-exponential-cutoff (6 = 1, solid curve) shape is a good fit 
for PSR J1057-5226, it does not fit the high energy points properly for PSR J1709-4429 
or PSR J1952+3252. A gradual-exponential-cutoff (6 < 1, dashed curve) model for these 
two pulsars agrees very well with all the points up to highest energies, consistent with 
the preference for this model over a 6 = 1 model in the full-band fit. A best-fit phase 
averaged s pectral shape with 6 < 1 is similar to what has been found for the Vela pulsar 



spectrum (lAbdo et all l2010fl ). Although a power law with a gradual cutoff, 6 < 1, does 
not have any physical meaning, it is a possible indication for large variations in -E C utofr with 
phase. The phase-averaged spectrum can be a combination of several spectra exhibiting 



simple exponential cutoff (6 = 1) with different cutoff en ergies. The same behavior was 
demonstrated in Vela using simulations (lAbdo et al.ll2010fl ). 



4-2.2. Phase-resolved Spectrum 

The large number of photons collected by the Fermi LAT from these pulsars allows us 
to obtain phase-resolved spectra in order to study the variation of the spectral parameters 
with phase. The phase bins for the phase-resolved spectra were defined such that in each bin 
a fixed-count of 500 photons for PSR J1057-5226, 1500 photons for PSR J1709-4429, and 
1000 photons for PSR J1952+3252 were collected in the energy- dependent radius defined in 
Section 14.11 The spectral analysis data set covering the 10° ROI around each pulsar was 
then divided into these phase bins, and an independent spectral fit was performed for each 
phase bin. For each spectral fit, the diffuse background and nearby sources were modeled 
with the same best-fit model obtained from the off-pulse analysis, thus only variations from 
bin to bin were due to the pulsar. The spectral fit parameters were obtained in each phase 
bin where the significance of the pulsar emission was above the 5a level. 

Although the data in each phase bin were not adequate for any of the pulsars to make a 
fit to test the shape of the cutoff while allowing b to vary, additional fits for PSR J1709— 4429, 
the brightest of these three pulsars, were performed in each phase bin fixing b to four different 
values between 0.5 and 1.0, and it was found that none of those gradual-cutoff-shape models 
gave a significantly better fit in any phase bin than the model with simple-exponential-cutoff 
shape {b = 1), with the current statistics, as opposed to the case of the phase-averaged 
spectra. The comparison of the cutoff energies obtained for each fixed b value also revealed 
that the .EcutofT and b parameters in Equation flTJ are strongly correlated (see Figur dTUj) . as 
also evident in the phase-averaged spectrum results with b fixed at 1 or left free (Table EJ. 
Thus, for these limited-counting-statistics spectra where we can not obtain a significantly 
better fit in any phase bin among different values of b, we assumed a shape of a power law 
with simple-exponential cutoff (b = 1) for each pulsar, allowing only the K, T, and E cuto Q 
parameters to vary in the fit. 

Figures [8] [10] and [12] show the variation of the spectral parameters with phase for each 
pulsar, while Tables [T] — [9] report the spectral parameters obtained from the fit in each phase 
bin for each pulsar. The quoted fluxes are normalized by the width of the phase bin. 

The spectrum of PSR J1057— 5226 was measured in eight bins in the phase range <p = 
0.28 — 0.64. No strong correlations between spectral parameters and phase are observed in 
this region. More data will be needed to fully characterize any dependence of the spectral 
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parameters with phase. 

The spectrum of PSR J1709— 4429 was measured in 15 bins in the phase range (p = 
0.085 — 0.624. The photon index varies slowly throughout the pulse profile, with a minimum 
around the leading edge of P2 and increasing toward the tails of the pulsed emission region. 
The cutoff energy shows a more drastic evolution, changing from a minimum of 1.14 ± 
0.23 GeV at the leading edge of PI to a maximum of 5.74 ± 1.04 GeV around the location of 
P3 seen in the pulse profile above 3 GeV. The cutoff energy has an increasing trend from PI 
to P2, but has a sharp maximum near P3, and falls off at the trailing edge of P2. 

Lastly, the spectrum of PSR J1952+3252 was measured in five bins at the locations of 
the peaks. The photon index is quite constant over the two peaks, while the cutoff energy 
has an increasing trend from PI to P2. 

Although more data are necessary to clearly assess any trends, the evolution of cutoff 
energy with phase for these pulsars seems to be compatible with the prediction that the 
b < 1 preference in the phase-averaged spectrum results from large variations of cutoff 
energy with phase, while b consistent with 1 implies that the cutoff energy should be constant 
with phase. The cutoff energy for PSR J1057— 5226 is quite constant over the phase, and 
the 6 = 1 cutoff shape is a sufficiently good fit for the phase-averaged spectrum. On the 
other hand, PSR J1709— 4429 shows the most drastic evolution of the cutoff energy over its 
phase and the best-fit phase-averaged spectrum of this pulsar requires the lowest value of 
b = 0.506 ± 0.021 stat ± 0.035 svs . 



4.3. Search for Off-pulse Emission 

The off-pulse phases of emission from these pulsars were also analyzed to search for 
potential unresolved nebular sources around the pulsars. For these searches, a minimum 
emission phase range in the off-pulse regions was selected to avoid contamination from the 
pulsed emission peaks. Spectral analyses, as described in Section 14.2} were performed in 
these off-pulse phase ranges by first assuming no steady emission exists at the location of 
the pulsar (null hypothesis) and then assuming a point source with a power-law emission 
spectrum at the pulsar locations (test hypothesis). The log-likelihood values of both fits for 
each pulsar are compared to calculate the test statistic (TS) of a possible off-pulse source 
at the location of the pulsar from the difference between the log likelihood with (Li) and 
without (L ) the source, i.e., TS = 2(Li — L ). 

We selected photons from the off-pulse region (<£> = 0.0 — 0.25) and ((p = 0.65 — 1.0) 
to search for nebular emission from PSR J1057— 5226. From the spectral fit performed with 
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gtlike, we cannot claim any clear detection of a nebula with a TS value of 15. We set an upper 
limit at 95% confidence level for the flux from the nebula to F(> lOOMeV) = 6.67 x 10~ 8 
ph cm _2 s~ 1 . 

For PSR J1709— 4429, we selected an off-pulse phase region of (if = 0.0 — 0.1) and 
(v? = 0.7 — 1.0) to search for a steady emission. We did not detect any significant point-like 
emission in that off-pulse phase interval with a TS of 23.4 for a test source at the location 
of this pulsar. We set an upper limit at 95% confidence level for the flux from the nebula of 
F(> lOOMeV) = 8.51 x 10~ 8 ph cm~ 2 s~ 1 . This result is consist ent with the non-detection of 



a poi nt-like off-pulse source by EGRET ( [Thompson et al.lll996l ) and HESS ( lAharonian et al. 



20051 ) around this pulsar. The analysis to search for extended emission around this pulsar is 
ongoing using the LAT data in the light of detection of a significant extended emission by 
HESS at an offset of 0.5° from the pulsar. 

The off-pulse emission from PSR J1952+3252 has been investigated in the phase interval 
of if = 0.75 — 1.00. We did not detect any significant point-like emission with a TS of 2.8 for 
a test source at the location of this pulsar. Thus, no emission is detected from CTB80, and 
we calculated an upper limit of F(> lOOMeV) = 6.5 x 10~ 8 ph cm _2 s _1 at 95% confidence 
level. 



5. Discussion 



The superior sensitivity and high statistics obtained with the LAT data allow con- 
struction of light curves and spectra with better precision, providing tighter observational 
constraints on the pulsed emission models than previou sly possible. Currently, there are two 
classes of high-energy pulsed emission models: the PC ( IDaugherty fc Hardinglll996l ) models 
predicting a lower altitude emission f rom near the magnetic poles of the neut ron star and 
high- altitude emission models like OG (Romani|[l996) or two-p ole caustic (TPC) ( Dyks & Rudak 
20031 ) models, like slot gap (SG) (jMuslimov &: Hardingll2004l ). which predict emission higher 
in the magnetosphere extending up to the light cylinder. 

The on-pulse regions of the light curves cover a wide range in phase, from 30% in the 
case of PSR J1952+3252 to 46% for PSR J1709-4429, while the radio emission peaks in a 
very narrow angle for all of them. This feature suggests that the gamma-ray beams cover a 
large solid angle compared to the radio emission and favors the prediction of high-altitude 
emission models, such as TPC and OG. 

The maximum energy e max of the pulsations provides a lower limit to the altitude of 
gamma-ray emission, since it must lie below the 7 — B magnetic pair creation absorption 
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threshold. This bound can be simply expressed in terms of the surface polar magnetic field 
strength lO^ -B^G, the pulsar period P (in seconds), and e max , inverting Equation (1) of 
BarineJ (120041 ): 

r > (e max 5i 2 /l.T6GeV) 2 / 7 P' 1 ' 1 R, (2) 

where i?* is the radius of the neutron star. The good photon statistics permit tracing the 
spectrum deep into the turnover for each of these pulsars: a super-exponential 7 — B pair 
attenuation feature is excluded by the spectra at energies above the cutoff energy E cntoS 
specified in Equation (JTJ. Accordingly, it is safe to adopt a value of around 2E cuto s for the 
maximum energy e max in bounding the emission altitude. If one chooses e max = 3.0 GeV in 
PSR J1057-5226 (Cutoff ~ 1-50 GeV), this yields r > 1.5R*. For PSR J1709-4429, the 
higher energy choice of e max = 8.9 GeV yields the bound r > 3.0-R*. For PSR J1952+3252, 
the limit is r > 1.8/2* when adopting e max = 5.6 GeV for -E cuto ff ~ 2.8 GeV. Clearly these 
bounds exclude emission very near the stellar surface, adding support for an SG or OG 
acceleration locale in these pulsars. 

The gamma-ray efficiency, r? 7 , of the pulsar can be obtained from the relation, r] 1 = 
L^/E, where L 1 is the total gamma-ray luminosity given as 



L 7 = 477/0(0;, (e) d 2 G 



100- 



(3) 



In this equation, Gi o is the observed energy flux above 100 MeV at the Earth line of sight (at 
angle (e to the rotation axis), d is the pulsar distance, and fa is the be aming correction factor 
that depends on the geometry of the emission. For the PC model (iDaugherty fc Harding 
19961 ). the gamma-ray emission originates at a few stellar radii from the surface, implying 
an emission with a small solid angle, that is, fa <C 1. For both OG and SG models, where 
the emission is far away from the neutron star, the resulting fa values can be near or even 
greater than 1. Aside from the unknown factor fa, a large uncertainty on the gamma-ray 
efficiency arises from the determination of the distance. The distances obtained from the 
literature for the three pulsars under study are listed in Table [1] and the observed energy 
flux values, G100, from these pulsars are listed in Table |6j 

Using these parameters in Equation (j3J), we found L 7 = (1.7 ±0.7) x 10 34 faergs -1 and 
rj y = (0.58 ±0.23) fa for PSR J1057-5226. A similar calculation for PSR J1952±3252 leads 
to L 7 = (6.9 ± 2.9) x lO^/nergs" 1 and an efficiency of r/ 7 = (0.019 ± 0.008) fa- Finally, 
for PSR J1709-4429 the L 7 is found in the range (2.58 x 10 35 fa -1.70 x 10 m fa) given 
the range of its estimated distances and its gamma-ray efficiency is in the range (0.076 fa — 
0-50 fa). From the EGRE T observations, there was a noted trend that gamma-ray pulsar 
effi ciency scales a s E 1 / 2 ( Thompson et all Il999l ). with a physical origin earlier predicted 
by iHardingl (Il98ll ). Fermi puls ars have confirmed this trend, at least for E > 10 34 , that 
?7 7 ~ 0.034( J E/10 36 ergs" 1 )- 1 / 2 (Ubdo et al.ll2010eh . PSR J1952±3252 follows this relation 
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almost exactly, while the efficiency of PSR J1057— 5226 lies above by a factor of 2, and the 
efficiency range of PSR J1709— 4429 lies substantially above the efficiency predicted with 
this relation. 

The number of main peaks and the peak separation may be used to constrain the 
geometry of the emission for high-altitude models. The measurements of the viewing angle, 
C, and the magnetic inclination angle, a, obtained from the torus fits to the X-ray images 
(INg fc Romanill2008l ) or from polarization measurements (jWeltevrede fc Wrightll2009l ). listed 
in Table [TU1 when a vailab le, provide additional constraints on the model, as described in 
Romani &: Watterj (120101 ). The possible ranges of geometrical parameters, a and (, can 
be obtained from the allowed phase space of the predicted geometry of the two geometrical 
models constrained with the measurements. The light curves modeled for this geometry are 
then compared with the observed light curves to determine a goodness of fit, \. The inferred 
angles, a and (, and the goodness of the fits , x> f° r the TPC and OG models, obtained 
from this procedure (IRomani fc Wattera |2010| ). are listed in Table [TT1 These geometrical 
considerations might give hints about the underlying emission model. 



Romani fc Watterd (120101 ) find that, in some cases, several models can produce accept- 



able fits. PSR J1952+3252 has only limited X-ray/radio constraints and both TPC models 
and OG models can produce light curves with comparable fit statistics. In the case of PSR 
J1709-4429, fits consistent with the externally measured a and £ can be found for both TPC 
and OG, although the fit statistic for the latter is substantially better. For PSR J1057-5226 
no viable TPC models are present anywhere near the radio-determined a and (, while the 
OG produces reasonable light curves for these values. 

The predictions for the pulse d emission geometry in these models assume the efficiency 
relation (w =)r] ~ (10 33 /£) 5 of lWatters et all (120091 ). and the gamma-ray efficiencies (or 
gap widths) derived from this assumption are also listed in Table [TTJ Combining these 
assumed gamma-ray efficiency values with the gamma-ray efficiencies derived from the ob- 
served gamma-ray luminosities, one can derive a flux correction factor fa for each pulsar and 
compare this value against the predicted flux correction factor for each pulsar from these 
models. For PSR J1952+3252 the flux correction factor fn, derived from this relation, is 
compatible with the predictions of both models. However, the derived flux correction fac- 
tor f n for the other two pulsars, PSR J1057-5226 and PSR J1709-4429, are significantly 
lower than the predictions of the geometrical models. Thus, those predictions of geometrical 
models for flux correction factors must be revised. In particular, PSR J1709— 4429 does 
appear to have a substantially higher efficiency than expected for its estimated distance and 
predicted beaming. 



The phase lag between the radio pulse and the first gamma-ray pulse is listed in the 
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third column of Table [10J These are larger than expected from a basic geomet ric model of 
both OG and two-pole models like the ones discussed in IWatters et al.l ( 120091 ). The main 
reason for the difference is their underlying assumption that the radio emission comes from 
the PC region near the surface. Ho wever, if an additional assumptio n of high-altitude radi o 
emi ssion is invoked as d iscussed in iKarastergiou fc Johnston! (120071 ) , IWatters et al.l ( 120091 ) , 
and IVenter et al.l ( 120091 ). the radio pulse will arrive at earlier phase because of aberration, 
and retardation due to the finite-time-of- flight effect, and the larger resulting phase lag with 
the gamma-ray pulse better matches the observed ones. 

The phase-resolved spectral measurements that are now achievable with the LAT will 
provide another important diagnostic for model comparisons. Among the three pulsars 
analyzed in this paper, only PSR J1709— 4429 so far has measurements with small enough 
errors and in fine enough phase bins to show significant variations. The variations seen for 
PSR J1709— 4429 i n photon index and cutoff energy s how patterns that a re similar to those 
measured for Vela ( lAbdo et al.l 1201 Off ) and Geminga ( jAbdo et al.l l2010dl ). notably a slight 
hardening in index in the bridge region and a rise to maximum values of the cutoff energy in 
the peaks. In general, the stronger variation of cutoff energy with phase may well be tied to 
variations in emission radius and/or magnetic-field-line curvature since most models predict 
that the emission at GeV energies is due to radiation-reaction-limited curvature radiation 
with exponential cutoffs in the range of the measured cutoff energies. 



6. Conclusion 



A common thread that links PSR J1057-5226, PSR J1709-4429, and PSR J1952+3252 
is that they were all detected by EGRET, although not nearly as strongly as Vela, Crab, and 
Geminga. One question is why thes e three were visible to EGRET when we now know there 
are at least 46 gamma-ray pulsars (lAbdo et al.ll2010el ). a number of which are in EGRET 
unidentified source error boxes. These three were not among the brightest Galactic sources 
seen by EGRET; in fact, PSR J1952+3252 did not appear in the third EGRET catalog, 
being detected only by its pulsations. The basic answer appears to be that these were 
among the brightest of the gamma-ray pulsars that had good radio timing in the EGRET 
era. Three other pulsars that had good radio timing were seen by EGRET as margina l 
detections (now all confirmed by LAT obser vations) : PSR J02 18 +4232 dKuiper et al.ll2000r) . 
PSR B0656+14 faamanamurthv et al.fl99~6T ) and PSR B1046 -58 JKaspi et al.lboooh . PSR B0740-2£ 



had good timing in the EGRET era, but the LAT detection ( lAbdo et al 



201(3) shows that 



the gamma-ray flux is well below the EGRET sensitivity. Many of the radio pulsars that 
have been detected by Fermi had timing programs started only after the primary EGRET 



-21 - 



mission was over. Therefore, it seems clear that PSR J1057— 5226, PSR J1709— 4429, and 
PSR J1952+3252 are not fundamentally different from the other LAT-detected pulsars. 

The other key question addressed by this paper is whether the unique features suggested 
by the EGRET observations are supported by the LAT observations. The answer is mixed. 

PSR J1057— 5226 in the EGRET era appeared to have a very high efficiency, high enough 
to argue against models with large beams as being physically unrealistic. Improved distance 
estimates now put this pulsar closer to Earth, but its efficiency still appears to be higher 
than that of most gamma-ray pulsars. 

PSR J1709— 4429 was unique among the EGRET pulsars in having an energy spectrum 
well fit by a broken power law, rather than a single power law or a cutoff at high energy. 
This fit was driven by the indications of significant pulsed emission above 10 GeV, a feature 
that is confirmed by the LAT data. Although the exponential cutoff spectral form fits the 
LAT data, it requires a high-energy cutoff and it was not dramatically better than a broken 
power-law fit until more than 9 months of data were included in the analysis. However, the 
LAT best fit to the phase-averaged spectrum yields b < 1 so that the spectrum falls off more 
slowly than a pure exponential, as may be expected from a cutoff whose energy varies with 
pulsar phase. This shape may have allowed an EGRET fit with a broken power law. 

PSR J1952+3252 showed no evidence of a spectral cutoff in the EGRET data, extending 
out to the 10—30 GeV energy band. The LAT results show a clear exponential cutoff at 
a much lower energy, so this feature from the EGRET era is not confirmed. As noted, 
however, this pulsar had very limited statistics with EGRET (only two photons above 10 
GeV). PSR J1952+3252 does not stand out among the LAT pulsars in any respect. 

PSR J1057-5226 and PSR J1709-4429 are among the minority of LAT pulsars that 
do not show widely spaced double pulses. Several pulsars in the first LAT catalog have 
closel y spaced double pulses similar to these , notably PSR J0 007+7303 (the CTA1 pul- 



sar) flAbdo et al]l2008h and PSR J1509-5850 flAbdo et alJl2010ah . So although these shapes 



are not the most common, they should be predicted by successful outer magnetosphere mod- 
els. Continued study will be needed in order to derive detailed information about the pulsar 
magnetospheres for these two. 
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Pulsar 


BName 


P 


P 


E 


Age (a) 


-^surface 


Distance 






(ms) 


(io- 15 ) 


10 34 ergss -1 


(10 3 yr) 


(10 12 G) 


(kpc) 


PSR J1057-5226 


B1055-52 


197 


5.83 


3.01 


535 


1.09 


0.72 ± 0.2^ 


PSR J1709-4429 


B 1706 -44 


102 


93.0 


341 


17.5 


3.12 


1.4 - 3.6^ 


PSR J1952+3252 


B1951+32 


39.5 


5.84 


374 


107 


0.486 


2.0 ±0.6^ 



Table 1: Various Cha racteristic Parameters for the Three Pulsars. Values are taken from 
the ATNF database ([Manchester et al.l 120051 ) except the distance estimates, for which the 
references are listed s eparately. We ad opt the same distance estimates used for the first 
Fermi pulsar catalog (lAbdo et al.ll2010el ). 



a Characteristic age: P/2P 

b Estimated using dispersion measure (jCordes fc Lazioll2002i ; IWeltevrede fc Wrighdl2009l ). 



Koribalski et al.1 (|1995j ): iMcGowan et al.l (|2004l ) 



^Evaluated from the kinematic model ( Greidanus fc Strom 199C)I) 



Pulsar 


Observer 


Validity Period (MJD) 


TOA 


rms 


PSR J1057-5226 * 


Parkes 


54220.2-55071.9 


47 


62 fis 


PSR J1709-4429 


Parkes 


54220.6-55072.3 


49 


1255 //s 


PSR J1709-4429 * 


Fermi LAT 


54647.4-55074.6 


23 


378 fjs 


PSR J1952+3252 * 


Nangay 


54607.1-55075.9 


17 


228 fis 



Table 2: Timing Observations. The number of TOAs used to derive the timing solutions and 
the post-fit root-mean-square (rms) of the timing solution are also given. The (*) denotes 
the observations that give the timing solution used in this analysis. 
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Energy Band 


PI Location 


PI HWHM 


P2 Location 


P2 HWHM 


0.1 - 100 GeV 
>3.0 GeV 
1.0 - 3.0 GeV 
0.3 - 1.0 GeV 
0.1 - 0.3 GeV 


0.31±0.01 
0.32±0.02 
0.32±0.01 
0.31±0.01 
0.31±0.01 


0.04±0.01 
0.04±0.03 
0.03±0.01 
0.03±0.01 
0.05±0.01 


0.59±0.01 
0.59±0.02 
0.59±0.01 
0.59±0.01 
0.58±0.02 


0.03±0.01 
0.04±0.02 
0.02±0.01 
0.03±0.01 
0.04±0.02 



Table 3: The location and Gaussian half-width-at-half-maximum (HWHM) values for the 
leading and trailing peaks (PI and P2) of the PSR J1057— 5226 pulse profiles in every energy 
band. These parameter values were obtained by fitting the peaks with Gaussian functions 
over a constant background. 



Energ 


;y Band 


PI Location 


PI HWHM (Outer) 


PI HWHM (Inner) 


0.1 - 


100 GeV 


0.242±0.002 


0.027±0.002 


0.072±0.008 


>10.0 GeV 








3.0 - 


10.0 GeV 


0.240±0.004 


0.008±0.003 


0.101±0.044 


1.0 - 


3.0 GeV 


0.243±0.003 


0.020±0.002 


0.103±0.016 


0.3 - 


1.0 GeV 


0.243±0.003 


0.028±0.003 


0.072±0.016 


0.1 - 


0.3 GeV 


0.238±0.004 


0.029±0.005 


0.061±0.012 


Energ 


;y Band 


P2 Location 


P2 HWHM (Inner) 


P2 HWHM (Outer) 


0.1 - 


100 GeV 


0.492±0.003 


0.118±0.013 


0.032±0.003 


>10.0 GeV 


0.501±0.009 


0.037±0.013 


0.009±0.010 


3.0 - 


10.0 GeV 


0.489±0.005 


0.126±0.032 


0.017±0.003 


1.0 - 


3.0 GeV 


0.496±0.003 


0.140±0.066 


0.023±0.002 


0.3 - 


1.0 GeV 


0.488±0.005 


0.115±0.018 


0.037±0.004 


0.1 - 


0.3 GeV 


0.496±0.007 


0.098±0.030 


0.032±0.006 



Table 4: The location and Lorentzian half-width-at-half-maximum (HWHM) values for the 
first and second peaks (PI and P2) of the PSR J1709— 4429 pulse profiles in every energy 
band. These parameter values were obtained fitting the peaks with asymmetric Lorentzian 
functions over a constant background. 
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Energy Band 


PI Location 


PI FWHM 


P2 Location 


P2 HWHM (Inner) 


P2 HWHM (Outer) 


0.1 - 100 GeV 
>3.0 GeV 
1.0 - 3.0 GeV 
0.3 - 1.0 GeV 
0.1 - 0.3 GeV 


0.154±0.001 
0.156±0.008 
0.155±0.002 
0.155±0.002 
0.152±0.004 


0.086±0.005 
0.047±0.022 
0.034±0.005 
0.073±0.006 
0.127±0.016 


0.639±0.004 
0.650±0.004 
0.635±0.005 
0.640±0.004 
0.650±0.008 


0.119±0.011 
0.097±0.017 
0.122±0.012 
0.112±0.015 
0.126±0.027 


0.034±0.007 
0.006±0.002 
0.029±0.007 
0.035±0.010 
0.017±0.007 



Table 5: The location and Lorentzian full-width-at-half-maximum (FWHM) and half-width- 
at-half- maximum (HWHM) values for the first and second peaks (PI and P2) of the 
PSR J1952+3252 pulse profiles in every energy band. These parameter values were obtained 
fitting the peaks with asymmetric Lorentzian functions over a constant background. 



PSR 


Prefactor (K) 
(1(T 10 phcm^s^MeV" 1 ) 


r 


-^Cutoff 

(GeV) 


b 


^100 

(10~ 7 phcm- 2 s _1 ) 


G*ioo 

(10" 4 MeVcm" 2 s" 1 ) 


J1057-5226 


1.25 ±0.07 ±0.05 


1.20 ±0.05 ±0.11 


1.50 ± 0.09±g;£jj 


1 (fixed) 


3.36 ±0.10 ±0.20 


1.77 ±0.01 ±0.16 


J1057-5226 


1.26 ±0.03 ±0.30 


1.19 ±0.11 ±0.32 


1.48 ±0.441°;^ 


0.99 ±0.15 ±0.07 


3.36±0.11±0.25 


1.77 ±0.01 ±0.16 


J1709-4429 


2.64 ±0.05 ±0.15 


1.71 ±0.02 ±0.06 


4.45 ± 0.23^8? 


1 (fixed) 


15.97 ± 0.24 ± 1.12 


6.98 ±0.08 ±0.90 


J1709-4429 


8.19 ±0.96 ±0.46 


1.30 ±0.04 ±0.06 


0.57 iO.lllg;^ 


0.51 ±0.02 ±0.04 


15.23 ±0.24 ±0.99 


6.99 ±0.09 ±0.71 


J1952+3252 


0.39 ±0.03 ±0.02 


1.57 ±0.07 ±0.15 


2.80 ± 0.37±g;|| 


1 (fixed) 


1.89 ±0.12 ±0.10 


0.92 ±0.03 ±0.08 


J1952+3252 


1.24 ±0.37 ±0.08 


1.13 ±0.31 ±0.12 


0.50 ± 0.02±g;iS 


0.57 ±0.03 ±0.06 


1.78 ±0.15 ±0.08 


0.91 ±0.16 ±0.09 



Table 6 Spectral parameters obtained from fitting the energy spectra with a power-law with exponential cutoff. The 
first four parameters are as defined in Equation flTJ and the last two parameters are the integrated photon flux and the 
integrated energy flux above 100 MeV. The first errors are statistical and the second are systematic errors calculated 
from the bracketing IRFs. 
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in 

Y^min 


in 

ymax 


±111.1x3 A. 


Z7> 

-C'cutoff 

l^ev ) 


f(~> i nn MpVI ( a ) 

1 \^ 1UU IVlC V J 


0.28 


0.32 


1.26 ± 0.11 


1.21 ± 0.19 


13.56 ± 1.13 


0.32 


0.37 


0.99 ± 0.12 


1.12 ± 0.15 


10.04 ± 0.88 


0.37 


0.41 


1.17 ± 0.12 


1.42 ± 0.22 


10.61 ± 0.96 


0.41 


0.45 


1.02 ± 0.11 


1.33 ± 0.18 


10.96 ± 0.95 


0.45 


0.49 


0.99 ± 0.12 


1.28 ± 0.19 


11.80 ± 1.05 


0.49 


0.53 


1.31 ± 0.10 


1.71 ± 0.28 


12.77 ± 1.08 


0.53 


0.58 


0.91 ± 0.12 


1.29 ± 0.17 


9.01 ± 0.82 


0.58 


0.64 


1.26 ± 0.11 


2.04 ± 0.37 


7.26 ± 0.71 



Table 7: Phase-resolved spectral parameters for PSR J1057— 5226. Only the statistical errors 
are quoted. 



"Normalized to the full phase range. 



Vmin 


Vmax 




Index 


-^cutoff 


F(> 100 Me V) ^ 










(GeV) 


(10- 7 cm- 2 s" 1 ) 


0.085 


0.165 


2 


24 ± 0.42 


1.68 ±1.58 


3.49 ±0.91 


0.165 


0.216 


1 


57 ±0.15 


1.14 ±0.23 


15.25 ± 1.21 


0.216 


0.242 


1 


73 ± 0.07 


2.94 ±0.50 


47.16 ±2.14 


0.242 


0.265 


1 


57 ±0.06 


3.42 ±0.48 


52.43 ±2.23 


0.265 


0.292 


1 


58 ± 0.06 


3.71 ±0.56 


41.56 ± 1.86 


0.292 


0.324 


1 


70 ± 0.06 


5.74 ± 1.04 


36.06 ±1.69 


0.324 


0.354 


1 


63 ± 0.06 


5.48 ±0.94 


33.99 ±1.66 


0.354 


0.385 


1 


56 ± 0.06 


3.86 ±0.58 


33.01 ± 1.57 


0.385 


0.416 


1 


48 ± 0.06 


3.69 ±0.53 


32.83 ±1.55 


0.416 


0.442 


1 


55 ± 0.06 


3.87 ±0.57 


43.89 ± 1.92 


0.442 


0.466 


1 


56 ± 0.06 


4.23 ±0.60 


46.22 ± 1.98 


0.466 


0.489 


1 


61 ±0.05 


4.74 ±0.71 


54.43 ± 2.22 


0.489 


0.514 


1 


65 ± 0.05 


4.71 ±0.73 


46.56 ± 1.99 


0.514 


0.554 


1 


75 ± 0.08 


3.18 ±0.67 


25.90 ± 1.49 


0.554 


0.624 


1 


72 ± 0.23 


1.34 ±0.51 


7.15 ±0.90 



Table 8: Phase- resolved spectral parameters for PSR J 1709— 4429. Only the statistical errors 
are quoted. 



a Normalizcd to the full phase range. 
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V-^min 


V-'max 


Index 


-^cutoff 

(GeV) 


F{> 100 Me V) ^ 
(lO^cm" 2 s" 1 ) 


0.094 


0.155 


1.69±0.15 


2.34±0.76 


5.57±2.30 


0.155 


0.224 


2.63±0.15 


1.67±0.44 


5.65±0.58 


0.490 


0.559 


1.32±0.21 


2.80±0.93 


2.61±0.43 


0.559 


0.617 


1.57±0.11 


3.68±0.93 


6.72±0.52 


0.617 


0.673 


1.59±0.09 


3.84±0.84 


7.86±1.79 



Table 9: Phase-resolved spectral parameters for PSR J1952+3252. Only the statistical errors 
are quoted. 



"Normalized to the full phase range. 



Pulsar 


A a 


S a a(°) 


C(°) 


PSR J1057- 


-5226 


0.28 ±0.03 


0.31 ±0.03 75 b 


Q9 b 


PSR J1709- 


-4429 


0.250 ±0.004 


0.242 ± 0.002 


53 c 


PSR J1952- 


-3252 


0.485 ±0.017 


0.154 ±0.001 





Table 10: The measured geometric parameters for these pulsars. The parameters derived 
from the fits to the pulse profile, the peak separation and phase lag between the radio and 
the first gamma-ray pulse, are listed in the second and third columns. The a parameter, 
obtained from the radio RVM for PSR J1057— 5226, is given in the fourth column and the 
£ parameter, as measured from torus fits to the X-ray image for PSR J1709— 4429 and from 
the RVM fit for PSR J1057-5226, is given in the last column. 



"Derived from the fits on the pulsed profiles discussed in Section |4~T1 
b Measured from radio RVM fit (jWeltevrede fc Wrightl[2009l) 
c Measured from torus fits to X-ray images (|Ng fc Romanill2008T ) 
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Pulsar 


w a 


a T pc(°) 


Ctpc(°) 


x a 


/o,tpc 


«og(°) 


Cog(°) 


x a 


/n,OG 


PSR J1057-5226 


0.183 


75 


69 


414 


0.91 


75 


69 


106 


0.82 


PSR J1709-4429 


0.017 


36 


53 


155 


1.30 


36 


56 


30 


0.89 


PSR J1952+3252 


0.016 


71 


84 


19 


1.10 


66 


78 


14 


0.84 



Table 11: The pre dicted geometr i c para meters. The gap width w from the assumption 
of Equation (1) in IWatters et al.l ( 120091 ) is given in Column 2. The predicted geometric 



parameters, a and ( and fn, for the high altitude model s, two-pole caustic 



'PC) and outer 



gap (OG), and the goodness of fit parameter x are from lRomani &: Wattersl (120101 ) and they 
are listed in Columns 3 — 10. 



a \, the fit statistics as defined in 



Romani fc Wattersl (|20 10T) indicates the goodness of the fits. 
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Fig. 1. — Energy dependence of the PSR J1057— 5226 pulse profile. The top panel is the 
pulse profile in the full energy band; the following four panels show the pulse profile in five 
different energy bands: 0.1 - 0.3 GeV, 0.3 - 1.0 GeV, 1.0 - 3.0 GeV, and > 3.0 GeV. The 
bin widths are 0.01 in phase except the highest energy band which is 0.02. Two rotation 
cycles are shown, and the fitting functions (thick solid lines) are superimposed on the light 
curves in the first cycle. The bottom panel shows the pulse profile at radio wavelengths, in 
arbitrary units (au), at 1.4 GHz provided by Parkes radio telescope. 
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Fig. 2. — Evolution of leading peak to trailing peak ratio with energy for PSR J1057— 5226. 
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Fig. 3. — Energy dependence of the PSR J1709— 4429 pulse profile. The top panel is the pulse 
profile in the full energy band; the following five panels show the pulse profile in five different 
energy bands as indicated. The bin widths are 0.01 in phase except the highest energy band 
which is twice this width. Two rotation cycles are shown, and the fitting functions (thick 
solid lines) are superimposed on the light curves in the first cycle. The bottom panel shows 
the pulse profile at radio wavelengths at 1.4 GHz provided by the Parkes radio telescope. 
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Fig. A.— Evolution of the PI to P2 ratio with energy for PSR J1709-4429. 
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Fig. 5. — PSR J1952+3252 light curves. The top panel is the gamma-ray pulse profile 
for E >0.1 GeV within the energy-dependent 68% containment radius around the pulsar 
position. Each bin is 0.01 in phase. Two rotation cycles are shown and the fitting functions 
(thick solid lines) are superimposed on the light curves in the first cycle. Four following 
frames: light curves in the indicated energy ranges. The bottom panel is the radio profile at 
1.4 GHz obtained at the Nangay radio telescope. 
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Fig. 6.— Evolution of the PI to P2 ratio with energy for PSR J1952+3252. 
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Fig. 7. — Gamma-ray spectrum for PSR J1057— 5226. The curve represents the fit with 
simple-exponential cutoff(6 = 1) shape in the full energy range of 100 MeV— lOOGeV. The 
spectrum points on the curve were obtained from an independent fit in each energy bin, with 
a model-independent method, as explained in the text. 
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Fig. 8. — Evolution of the photon index (top panel) and the cutoff energy (bottom panel) of 
the pulsed emission spectrum through the pulse profile of the PSR J1057— 5226. 
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Fig. 9. — Gamma-ray spectrum for PSR J1709— 4429. The solid curve represents the fit 
with simple-exponential-cutoff {b = 1) shape and the dashed curve represents the fit with 
gradual-exponential-cutoff (b < 1) shape in the full energy range of 100 MeV- 100 GeV. The 
spectrum points on the curve were obtained from an independent fit in each energy bin, with 
a model- independent method, as explained in the text. 
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Fig. 10. — Evolution of the photon index (top panel) and the cutoff energy (bottom panel) 
of the pulsed emission spectrum through the pulse profile of the PSR J1709— 4429. 
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Fig. 11. — Gamma-ray spectrum for PSR J1952+3252. The solid curve represents the fit 
with simple-exponential-cutoff (6=1) shape and the dashed curve represents the fit with 
gradual-exponential-cutoff (b < 1) shape in the full energy range of 100 MeV- 100 GeV. The 
spectrum points on the curve were obtained from an independent fit in each energy bin, with 
a model-independent method, as explained in the text. 
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Fig. 12. — PSR J1952+3252: evolution of the photon index (top panel) and the cutoff energy 
(bottom) with phase. Only statistical errors are shown. 



